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Abstract. We analyze the properties of the quasiparticle excitations of metallic antiferromagnetic states in 
a strongly correlated electron system. The study is based on dynamical mean field theory (DMFT) for the 
infinite dimensional Hubbard model with antiferromagnetic symmetry breaking. Self-consistent solutions 
of the DMFT equations are calculated using the numerical renormalization group (NRG). The low energy 
behavior in these results is then analyzed in terms of renormalized quasiparticles. The parameters for 
these quasiparticles are calculated directly from the NRG derived self-energy, and also from the low energy 
fixed point of the effective impurity model. From these the quasiparticle weight and the effective mass are 
deduced. We show that the main low energy features of the fc-resolved spectral density can be understood 
in terms of the quasiparticle picture. We also find that Luttinger's theorem is satisfied for the total electron 
number in the doped antiferromagnetic state. 

PACS. 7 1.10.Fd,71.27.-Ha 



1 Introduction 

. The nature of the metallic antiferromagnetically ordered 

■ state in strongly correlated systems has been subject of 
] study for over two decades, but still remains to be fully un- 

■ derstood. Interest in this topic has been stimulated by the 
' fact that the high temperature superconductivity of the 
I cuprates emerges from the doping of an antiferromagnetic 

■ insulating compound, such as La2Cu04 [BIl]. The sim- 
[ plest models to describe the electrons in the Cu02 planes 
I of the cuprates are the two dimensional t- J- model and 
' Hubbard model. Much of the initial effort went into the 

study of a single hole state of these models in an antifer- 
' romagnetic background. For the motion of this hole, there 
, is a competition between the gain in kinetic energy from 
the hopping and its disruptive effect on the antiferromag- 
netic order, and consequent loss of potential energy. As a 
result a hole excitation becomes a quasiparticle or mag- 
netic polaron, heavily dressed by antiferromagnetic spin 
fluctuations (see review article by Dagotto |3j and refer- 
ences therein). Much of this work, however, relied on exact 
diagonalization or quantum Monte Carlo methods, which 
are limited to small clusters and very few hole excitations, 
and cannot be readily extended to study the many-hole, 
finite doping situation. 

More recent studies capable of describing finite dop- 
ing have concentrated on the relation between the antifer- 
romagnetic fluctuations and superconducting order (for 
a review see [1] and the references therein). One of the 
main motivations is to understand whether the exchange 
of these types of fluctuations can provide a purely elec- 



tronic mechanism for inducing superconductivity. Here, in 
this paper, we focus on the metallic antiferromagnetism, 
the doped state with long range antiferromagnetic order. 
Our interest is to examine how well the low energy ex- 
citations in this ordered state can be described in terms 
renormalized quasiparticles. To tackle this problem we use 
the infinite dimensional Hubbard model. 

The simplification in the infinite dimensional limit is 
that the electron self-energy becomes local in character, 
with no wavevector dependence . The self-energy then 
depends only on the frequency, as is the case for impurity 
models, allowing the lattice problem to be cast in the form 
of a self-consistent impurity model. There are several rea- 
sonably accurate techniques for solving this effective im- 
purity problem, a very accurate one for the zero and low 
temperature regime being the numerical renormalization 
group approach (NRG) [7]. 

Recently we studied the effect of a magnetic field on 
the quasiparticle excitations in the strong correlation regime 
of the infinite dimensional Hubbard model using the NRG 
method j8j . We also extended a form of renormalized per- 
turbation theory (RPT), originally developed for impu- 
rity models [9], to this model, and used it to calculate 
the local dynamic spin susceptibilities, obtaining results 
in good overall agreement with those from the NRG. In 
this paper we extend this combination of renormalization 
techniques, NRG and RPT within dynamical mean field 
theory, to look at the low energy excitations of the infi- 
nite dimensional Hubbard model in a staggered field, and 
in antiferromagnetic broken symmetry states. Extensive 
calculations of the antiferromagnetic states in the Hub- 
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bard model using the DMFT-NRG approach have already 
been reported in the paper of Zitzler, Pruschke and Bulla 
[To]. We confirm their results for the phase diagram and 
extend the calculations and analysis to the description of 
the low energy renormalized excitations, and how these 
can be described within the framework of a renormalized 
perturbation theory. 



2 Antiferromagnetic Broken Symmetry in 
DMFT 

In considering the response of the Hubbard model [11] to 
a staggered magnetic field and antiferromagnetic order, 
we take the case of a bipartite lattice, which consists of 
two sublattices A and B such that the nearest neighbors 
of a site in the A sublattice are on the B sublattice and 
vice versa. The Hamiltonian for the Hubbard model can 
be written in the form, 

i,j.(7 i,a 

^ ^(Mo-C^^j CT'^A,i,cr + t^-<y(^B,i,a^B, 1,(7)^ (1) 

i.cr 

where the hopping matrix element is taken as tij — —t 
between nearest sites i and j only, and zero otherwise, 
and a = A, B . A staggered field iJ* 



m = 



H for i E A sublattice 
—HioriEB sublattice 



(2) 



has been included so that ii^ = ^J^+crh, where h ~ g^^H /2 
with the Bohr magneton /^b- The non- interacting part 
of the Hamiltonian -ffo,M can be diagonalized in terms of 
Bloch states and then expressed in the form. 



(3) 



where C^ ^. = (c^.fc.a' '^L.fc.o-)' ^^^^ ^he matrix Mk^cr is 
given by 



(4) 



The k sums run over a reduced Brillouin zone, and the 
energy of the Bloch state is £k — J2j ^^j'^^^^''^^'' ■ The 
free Green's function matrix ^{uj) is given by (w — 

Mk^a-)~^- The poles of the free Green's function give the 
elementary single particle excitations, which are given by 



El^iU = 0) = -tioih)±^h^+el 



(5) 



where fJ.o{h) is the chemical potential of the noninteracting 
system in a staggered field. This illustrates that the elec- 
tronic excitations are split into two subbands for a finite 
staggered field. 

Notice that we have adopted a special choice of ba- 
si^ {'^A k tn'^B k a} ^^''^ [T2l[T0] . Another common basis to 



study antiferromagnetic and spin density wave symmetry 

}, where Qq is the recip- 



(7 ' k-^qojCT 



(SDW) breaking is {c^ 

rocal lattice vector for commensurate SDW ordering. The 
bases can be related by a linear transformation. 



'-k.a 
k+qa.a 



^'A,k,(y 
^B,k,cT 



(6) 



For the latter basis the matrix M^^^ would be diagonal 
in the kinetic energy term and the symmetry breaking 
is offdiagonal. For our study in the DMFT framework the 
A— i?-sublattice basis is, however, more convenient and we 
will use it throughout the rest of this paper. It is possible, 
of course, to relate the obtained quantities with the help 



of dll) to the {Cfc 



<TJ '^k+qo,<T 



} basis. 



We can generalize the equations to the interacting prob- 
lem by introducing a self-energy Ea^k^a-i^), so that the 
matrix Green's function can be written in the form 



1 



tB,k,<ri.^) 



CA,fc,cr('^)CB,fc,o-(w) 



-£k 



-£k C,A,k,a{^^)j 



(7) 

where Ca,fc,o-(w) As we are dealing 

with the infinite dimensional limit of the model, we take 
the self-energy to be local so we can drop the k index. This 
is the reason why the self-energy has a single site index 
a — A,B and no offdiagonal terms appear in equation fT]) . 
The symmetry of the bipartite lattice gives Es^ai^) — 
Z'^^_cr(w) = E-cri(^) and hence 

where we have simplified the notation. To determine these 
quantities Eaito) it is sufficient to focus on the A sublattice 
only. 

Summing the first component in the Green's function 
in equation ([7]) over k we obtain the Green's function for 
a site on the A sublattice, G^°'^{u!), 



de 



Po(£) 



(8) 



where po{£) is the density of states of the non-interacting 
system in the absence of the staggered field. 

In the DMFT this local Green's function, and the self- 
energy Ecy{uj), are identified with the corresponding quan- 
tities for an effective impurity model [12j . This implies 
that the Green's function t/o,CT(w) for the effective impu- 
rity in the absence of an interaction at the impurity site 
is given by 



g,l{co)^Gi°%uj)-' + E,{u). 



(9) 



We can take the form of this impurity model to correspond 
to an Anderson model [13 in a magnetic field. 



Ham = ^Ed^adld^ + Und^^nd,i 

a 



(10) 
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where Ed, a = Ed — crh is the energy of the locahzed level at 
an impurity site in a magnetic field H, U the interaction at 
this local site, and Vk^a the hybridization matrix element 
to a band of conduction electrons of spin a with energy 
Efe^cr- As we are focusing on an A site as the impurity we 
take H — Hs- 

The one-electron Green's function for the impurity site 
of this model is given by 



1 



where 



Kc,{ijj) ~ ' 



(11) 



(12) 



If this impurity Green's function is equated to the local 
lattice Green's function G^°'^{uj), we identify e^a — —fJ-a 
and from equation ([9]), /•^^(cj) is given by 



K^{uj) =uj + fi^ - Go,liuj) 



(13) 



The function Kcr{uj) plays the role of the effective medium 
and has to be calculated self-consistently. 

The self-consistent calculations for K^iuj) can usually 
be performed iteratively. Starting from a conjectured form 
for Ko-{llj), the NRG method is used to calculate the self- 
energy of the effective Anderson model, from which the 
impurity Green's function G™P(ti;) in (fTTjl and the local 
Green's function for the lattice G^°'^{uj) in ([5]) can be de- 
duced. If these two Green's functions do not agree, then 
equation ([9|) is used to derive a new starting value for 
Ka-{uj) and the process continued until self-consistency is 
achieved. 

To find antiferromagnetic solutions, we calculated self- 
consistent solutions for a decreasing sequence of staggered 
magnetic fields to see if broken symmetry solutions of this 
type exist as the staggered field is reduced to zero. For 
the non-interacting density of states pois) we take the 
Gaussian form po(e) = e^'^^* ^ /^t*, corresponding to 
an infinite dimensional hypercubic lattice. It is useful to 
define an effective bandwidth W = 2D for this density of 
states via D, the point at which po{D) = po(0)/e^, giving 
D = ^/2t* corresponding to the choice in reference [14]. 
In all the results we present here we take the value W = 
4. In the NRG calculations we have used the improved 
method |15pi6j of evaluating the response functions with 
the complete Anders-Schiller basis [17] , and also determine 
the self-energy from a higher order Green's function fTS]. 

In figure [1] we show the self-consistently calculated lo- 
cal spectral density for the spin-up (upper panel) and spin- 
down electrons (lower panel) at an A site with U — S and 
5% hole doping (from the state at half-filling) for various 
values of an applied staggered field. The staggered mag- 
netic field induces a sublattice magnetization. 



(14) 



so that these spectra are quite different. For this set of 
parameters, this difference persists as the staggered field 
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Fig. 1. (Color online) The spectral densities for the spin-up 
electrons (upper panel) and spin-down electrons (lower panel) 
at the A site for various values of the applied staggered field 
for (7 = 3 and x = 0.95 



is reduced to zero so that we have a spontaneous sub- 
lattice magnetization corresponding to spontaneous anti- 
ferromagnetic order. For the case away from half filling, 
(5 ^ 0, we have to keep adjusting the chemical potential 
when iterating for a self-consistent solution. It shows a 
slightly oscillatory behavior when iterating for a specific 
filling a;, and we follow the procedure described in refer- 
ence [TU] . This feature is related to the fact that the calcu- 
lations are for a metastable ground state and instabilities 
to more complicated ground states for antiferromagnetic 
ordering than the homogeneous, commensurate Neel state, 
which forms the basis for these DMFT calculations, can 
occur [T9l[20l[2Tl[22llMl[Mll25l[M][T0] . As far as phase sep- 
aration in the ground state is concerned, the results of 
our calculations are generally in line with the conclusions 
in |10j as they are carried out within the same frame- 
work. The focus of this work is, however, the analysis of 
generic quasiparticle properties in a doped antiferromag- 
netic state. We consider the approach as a valid, approxi- 
mate starting point for this endeavor, but modifications to 
the results presented here can occur for calculations based 
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on a more complicated ground states not accessible within 
the DMFT framework. For a more extensive discussion of 
the applicability of the DMFT in this situation we refer 
to the earlier work [10]. 

From results of this type of calculation, we have built 
up a global antiferromagnetic/paramagnetic phase dia- 
gram as a function of the doping S and the on-site inter- 
action U. This phase diagram is shown in figure [21 where 
the value of the corresponding sublattice magnetization is 
shown in a false color plot. We have added a line separat- 
ing the spontaneously ordered and paramagnetic regimes. 
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Fig. 2. (Color online) Phase diagram showing the doping and 
the I] dependence of the sublattice magnetization as de- 
duced from the DMFT-NRG calculations. 



At half filling (6 — Q axis) the spontaneous magnetiza- 
tion increases with V . We can see that the antiferromag- 
netic order from the half filled case persists when holes 
are added. The value of the critical doping 8c at which 
the antiferromagnetism disappears depends on the on-site 
interaction U . We expect that for small U the critical dop- 
ing (5c will increase with U since a tendency to order only 
appears when an on-site interaction is present. From the 
mapping to the t— J model we also expect that for large U 
the antiferromagnetic coupling J decreases and therefore 
the order is destroyed more easily. The values of U are, 
however, not large enough to display this trend. 

If we compare these results with the phase diagram 
given by Zitzler et al. jlOj we see that they are in very good 
agreement. In their case the antiferromagnetic region was 
mapped out to values oiU ~ 4.5. As the iterations tend 
to oscillate, as discussed before, there is a problem of ob- 
taining a self-consistent antiferromagnetic solution in the 
large V regime. We have managed to extend the diagram 
to somewhat larger values of U by stabilising the calcula- 
tions by averaging the effective medium over a number of 
iterations. 



3 Local Quasiparticle Parameters 

To examine the nature of the low energy excitations, we 
will assume that the self-energy is non-singular at 

w = so that, at least asymptotically, it can be expanded 
in powers of w. This assumption is not expected to be valid 
close to the quantum critical point when the magnetic 
order sets in, but to be a reasonable assumption otherwise. 
We also assume that the imaginary part of the self-energy 
vanishes which is confirmed by the numerical results of the 
DMFT-NRG calculations. We will retain terms to order 
Lo only for the moment. The higher order corrections will 
be considered later. We then find for C(t('^), 

CaH -co(i-i:;(o)) + Ai,-i:,(o) (15) 

= z^^(a) + /io,ff), (16) 

where 

li^,a = z„{[i-J:am. and - l-i;;(0). (17) 

The interacting Green's function ([7]) has poles at the roots 
of the quadratic equation. 



C,(^)C-.(co)-e| = 0. 
The solutions of this equation are 



All 



,-,2 



(18) 



(19) 



where = ^zf^efe, Ap. = (/io,T - Mo,i)/2, and /i = 
(Mo,t + Mo,i)/2. This has the same form as for the non- 
interacting system in a staggered field so we can in- 
terpret these excitations as quasiparticles coupled to an 
effective staggered magnetic field = Ajl/gp^, with jl 
playing the role of a quasiparticle chemical potential. This 
equation gives the dispersion relation for these single par- 
ticle excitations, which can be regarded as constituting a 
renormalized band, or bands as there are two branches. 
The term magnetic polaron is sometimes used to describe 
these single particle excitations in states of magnetic or- 
der, because of the analogy with the motion of a particle 
in a lattice to which it is strongly coupled, where the ex- 
citation is termed a polaron. 

The corresponding density of states of these free local 
quasiparticles on the sublattice is 



p, — aAfi 
fi -\- crA fi 



(20) 

for + //| > \Ajl\, and is zero otherwise. In the case of 
a half-filled band fi ~ Q and there is a gap at the Fermi 
level Ef = 0. 

To determine this local quasiparticle density of states 
in the presence of the symmetry breaking staggered mag- 
netic field we need to calculate and /io.o- for each spin 
type. Using the NRG we can do this in two ways. As the 
DMFT-NRG calculations give us the self-energy Saii-^) di- 
rectly, we only need its value, and that of its first derivative 
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at uj — 0, to deduce both and //q.o- using equation (|17p . 
However, because the model is solved using an effective 
impurity model, we can also deduce these quantities indi- 
rectly from the many-body energy levels of the impurity 
on approaching the low energy fixed point [57|. This sec- 
ond method gives us not only a check on the results of the 
direct method, but also allows to deduce some informa- 
tion about the quasiparticle interactions, as we shall show 
in the next section. 



3.1 Calculation of Renormalized Parameters 

To describe how the renormalized parameters are deduced 
from the energy levels of the NRG calculation, we need to 
outline how the NRG calculations are carried out. Fol- 
lowing the procedure introduced by Wilson [28] , the con- 
duction band is logarithmically discretized and the model 
then converted into the form of a one dimensional tight 
binding chain, coupled via an effective hybridization Va to 
the impurity at one end. In this representation K^{ll!) = 

\Va\'^ gQ^J {ll>) , where g^^J (uj) is the one-electron Green's 
function for the first site of the isolated conduction elec- 
tron chain of length N. The impurity Green's function for 
this discretized model then takes the form. 



to the poles of the quasiparticle Green's function given in 
equation (I24p . with values of Va and idcr, which are inde- 
pendent of iV as TV — + oo. We can test this hypothesis by 
substituting the values, w = Ep{N) and lu — E'^{N), into 
the equation. 



ed.-|Kp.g^;:\c.) = 0, 



(25) 



and deduce values of Va and Cdo-, which will in general 
depend upon N . From these we can deduce — Y^a/Va]"^ 
and /io.CT = — Edcr, which will also depend upon N , but if 
the lowest single particle excitations of the system do cor- 
respond to free quasiparticles, the values of Za and /io.o- 
will become independent of TV for large TV. It should be 
noted that we need both the particle and hole excitations 
for each spin to determine the four renormalized parame- 
ters. The parameters corresponding to spin up involve the 
particle excitations with spin up and the hole excitations 
with spin down. 

That parameters can be found, which are independent 
of TV for large TV, can be seen in figure[31 where we take the 
results of a K„{uj) and fx^ from the antiferromagnetic self- 
consistent solution for the Hubbard model with [/ = 3 and 
10% doping, using a value for the discretization parameter 
A = 1.8. 



(21) 



We can find the quasiparticle excitations of this model 
by expanding the self-energy Sa{u)) in the denominator of 
this equation to first order in a>, and write the result in 
the form. 



L0-eaa-WgZ\^) + O{u^^) 



(22) 



where 



ida = za[eda + S^m, |K|2 = z,|y,|2. (23) 



We can then define a free quasiparticle propagator, Go,o-(w), 
viz 

G'o7(^) = Z L.. iN). (24) 



(C) 



and interpret z„ as the local quasiparticle weight. 

In the NRG calculation the many-body excitations are 
calculated iteratively, starting at the impurity site, and 
increasing the chain length TV by one site with each itera- 
tion. When the matrices become too large to handle, only 
the lowest 500-1500 states are kept at each iteration. The 
many-body energy levels for the TVth iteration and the 
set of quantum numbers TVf, -Em (TV), depend on the chain 
length TV and the discretization parameter A > 1. When 
TV becomes large these energy levels go to zero as yl^^/^. 
We now conjecture that the lowest single particle E'^{N) 
and single hole excitations E^{N) determined from the 
NRG many-body excitations correspond to quasiparticle 
excitations. If this is the case then they should correspond 
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Fig. 3. (Color online) The TV-dependence of the renormalized 
parameters Za and /io,o- for U = 3 and x = 0.9. 



It can be seen that after about 25 iterations all the values 
deduced for z^ and /io,cr become independent of TV. In the 
next section, where we compare these results with the cor- 
responding values deduced directly from the self-energy, 
we get further confirmation that the values deduced re- 
ally do describe the quasiparticle excitations of the lattice 
model. 

When two or more quasiparticles are excited from the 
interacting ground state, there will be an interaction be- 
tween them. For the Anderson impurity model this inter- 
action will be local and can be expressed as f/, a renor- 
malized value of the original interaction of the 'bare' par- 
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tides. The value of U can be deduced by looking at low- 
est lying two-particle excitations derived from NRG cal- 
culation. These could either be two-particle excitations, 
EJ^^^{N), two-hole excitations, e\^^{N) or a particle- hole 
excitation E^lJ (N). By looking at the difference between 
a two-particle excitation and two single particle excita- 
tions, Ej,p{N) - El{N) - E^{N), as a function of N 
we can deduce an effective interaction U^^^{N) between 
these two quasiparticles, as has been described fully ear- 
lier for the standard Anderson model [27j . In a similar way 
we can deduce an effective interaction between two holes, 
^hh(^)' or a particle and hole, -Ul^i^)- To be able to 
define a single quasiparticle interaction U, not only must 
U^^^{N), UthW and C/Ja^(A^), give values which are in- 
dependent of N for large N, these values must be equal 



so c/;^^ 



1,1 

hh 



u. 




U, together with g^^{ijj), can be used as an alternative 
way of specifying this model. The renormalized perturba- 
tion theory (RPT) is set up by expanding the self-energy 
to order w, as earlier, but retaining all the higher order 
correction terms in a remainder term, 



Fig. 4. (Color online) The A''-dependence of the renormal- 
ized particle-particle, particle-hole and hole-hole interactions 
for [7 = 6 and x = 0.9, showing that they converge to a unique 
value U . 



In figure [4] we give the values of C/pV(^), ul'^{N) and 
Ulh{N) as deduced from DMFT-NRG calculation for the 
Hubbard model in an antiferromagnetic state with [/ = 6, 
10% doping and A = 1.8. It can be seen that the three 
sets of results settle down to a common value U . 

We can go further and identify U with the local quasi- 
particle 4- vertex interaction for the effective impurity model. 



[/ = ZT-^in,i.i,T(0, 0,0,0), 



(26) 



where /^^^^^^^^(a;!, iu;2, W3, 0-14) is the total 4-vertex at the 
impurity site, which is equal to the same quantity for a 
site in the lattice model. With this interpretation it is 
possible to identify these parameters with those used in a 
renormalized perturbation expansion. The parameters, V , 
£d,a and J7, together with ^{uj), specify the effective im- 
purity model. The renormalized parameters, V ^ ed,a and 



S,{u:) - ^,(0) + ujS'M + K'^'i^), 



(27) 



where i7^®™(w) is the remainder term. On substituting 
this into the equation for the impurity Green's function in 
equation pip , we can deduce a general expression for the 
quasiparticle Green's function in the form, 



G™P(^) 



1 



(28) 



where K„{lo) — z„K„{uj) and S^{u>) = ZaS]^™{ijj) plays 
the role of a renormalized self-energy. A diagrammatic 
perturbation theory can then be carried out for Sry{'jj) 
in terms of the free quasiparticle propagators, with ad- 
ditional diagrams arising from counter terms, which are 
required to prevent over-counting (renormalization condi- 
tions) pSlf^lSO] . This form of perturbation theory is valid 
for all energy scales but is particularly effective for calcu- 
lating the low energy terms arising from the quasiparticle 
interactions. For the symmetric Anderson impurity model 
it has been shown that this perturbation theory taken to 
second order in [/, gives the exact spin and charge suscep- 
tibilities at T = 0, and the exact contribution to the 
conductivity [9]. 

Because, within DMFT, the self-energy for the lat- 
tice is the same as that for the effective impurity, we can 
equally well use the effective impurity model to calculate 
it. This means that we can use the renormalized pertur- 
bation theory for the effective impurity model to estimate 
the correction terms to the free quasiparticle picture aris- 
ing from the quasiparticle interactions. 



3.2 Local Quasiparticle Weight 

We now consider the values of the local quasiparticle weight 
factor Zcr, commonly known also as the wavefunction renor- 
malization factor. This is an important factor in deter- 
mining the parameters needed to describe the low energy 
behavior of the system. When there is no fc-dependence 
of the self-energy, as is the case for infinite dimensional 
models and DMFT, the effective mass of the quasiparti- 
cles in the paramagnetic state is proportional to \/ We 
show later that in the antiferromagnetic state the expres- 
sion is more complicated and depends both on and the 
renormalized chemical potential /io,^- We have determined 
Za from the NRG results by the two methods described 
and give the values of z„ deduced for both spin types as a 
function of doping in figure [5l The results are for the case 
U — where there is antiferromagnetic order and the ex- 
ternal staggered field has been set to zero. It can be seen 
that there is a reasonable agreement between the values 
obtained by the two different methods of calculation. Visi- 
ble differences can be attributed to the inaccuracies when 
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Fig. 5. (Color online) The local quasiparticle weight z^, as de- 
duced directly from the self-energy and also from the impurity 
fixed point (FP) for [7 = 3 and various dopings. 



numerically computing the derivative of the self-energy, 
whose calculation involves a broadening procedure. When 
the system is doped but still ordered we have Z| ^ z^, and 
the renornialization effects are stronger for the minority 
(down) spin particles on the sublattice. This is similar to 
the results we found for a doped Hubbard model in a para- 
magnetic state in the presence of a strong uniform mag- 
netic field [5] . For a certain range of dopings the values of 
Z| and zi do not vary much. The tendency is that zj^ first 
decreases and later increases, whereas z-^ decreases over 
the whole range until both of them merge at the doping 
point where the antiferromagnetic order disappears. 

The results for the corresponding case with [/ = 6, a 
value which is larger than the bandwidth, are shown in 
figure [6l 
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Fig. 6. (Color online) The local quasiparticle weight Zcr as 
deduced directly from the self-energy and from the impurity 
fixed point (FP) for = 6 and various dopings. 



On the whole the behavior is quite similar to that for 
the case U — 3, only that the renormalization effects 



are more pronounced. For a range of dopings the local 
quasiparticle weights do not change much and have the 
same tendency as described above. The implications for 
the spectral quasiparticle weight and the effective mass 
enhancement will be discussed in detail later. 



3.3 Renormalized chemical potential 

In figure [7| we give the results for the renormalized chemi- 
cal potential, /lo.o- [defined in equation (fTT]) and ((23l) ]. for 
the two spin types in the spontaneously ordered antifer- 
romagnetic states for U = 3 and U = 6 for a range of 
dopings. 
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Fig. 7. (Color online) The renormalized chemical potential 
/io,CT as deduced directly from the self-energy and from the 
impurity fixed point (FP) for various dopings for [7 = 3 (upper 
panel) and U = 6 (lower panel). 



The values calculated by the two different methods can be 
seen to be in good agreement here, as well. We have added 
the values for the half filled case. These were calculated 
from the self-energy in the gap at cj = 0. The general 
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behavior of the values for /io^o- for the case with U — 6 is 
very similar to the case with smaller U 

The renormalized chemical potential /io.cr is an impor- 
tant parameter in specifying the form of the local sublat- 
tice quasiparticle spectral density From equation 

l(20|) it can be seen that, as a; — > — /2o,(t, Po,(t(w) behaves 
asymptotically as 

PoA^) - , \~ : (29) 

so the quasiparticle density of states has a square root 
singularity at w = —flo.a- On the other hand, however, as 
uj — /lo,-(T, Po,ct('^) behaves as 

/5o,<t(w) ^ \/^^ + /Uo,-CT, (30) 

so the quasiparticle density of states goes to zero at a; = 
^Ao.-cr- Between the two points, uj — —jlo.cr and lo = 
—flo^-cr, the quasiparticle density of states has a gap of 
magnitude 2Afl. As can be seen in figure[7]this free quasi- 
particle gap decreases with increasing doping and closes 
in the paramagnetic state. If we take into account the 
values at half filling we see a strong reduction of 2Ajl, 
when doping the system. We also see that /xo,| drops to 
small negative values for finite hole doping, which corre- 
sponds to the fact that the Fermi level then lies within the 
lower band. These features will be seen clearly in the fig- 
ures presented in the next section, where we compare the 
quasiparticle densities of states with the full local spectral 
densities calculated from the DMFT-NRG. 



3.4 The Quasiparticle Interaction 

The quasiparticles can be further characterized by an ef- 
fective interaction U as described before. In figure [8] we 
plot the doping dependence of the renormalized interac- 
tion over a range of dopings and U = S and U = 6. 
We can see that in both cases the values decrease with 
increasing doping. Hence, the effective quasiparticle inter- 
action is stronger for a smaller hole density. For a certain 
range of dopings, however, U does not vary much. We can 
also see that the ratio U/U for the effective interaction as- 
sume smaller values the larger the bare U becomes. Also 
the absolute value of U, i.e. without the scaling with U as 
in figure [H is smaller for larger bare U for the full range 
of dopings. This effect of smaller quasiparticle interactions 
for the stronger coupling case can be seen as sharper quasi- 
particle peaks for larger U as will be discussed in the next 
section. 



4 Spectra and Quasiparticle Bands 
4.1 Local Spectra 

In this section we examine how well the local sublattice 
quasiparticle density of states /5o,cr(w), evaluated from equa- 
tion (|20p with the renormalised parameters, describes the 
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Fig. 8. (Color online) The renormalised quasiparticle interac- 
tion U/U as deduced from the impurity fixed point for various 
dopings and U = 3,6. 



low energy features seen in the local spectral density p„ (w) 
calculated from the DMFT-NRG. At half filling there is 
a gap at the Fermi level, so there are no single particle 
excitations in the immediate neighbourhood of the Fermi 
level, and this is not a very interesting case to consider. 
We look in detail at the case of 10% doping where the 
Fermi level lies at the top of the lower band, and consider 
the two cases U = 3 and ?7 = 6. In the upper panel of 
figure [5] we compare the spectral density Pi{uj) with the 
corresponding quantity Z|Po,t('-^)i from the quasiparticle 
density of states. 

We see that the behavior near the Fermi level (w = 0), 
and the singular feature seen in the lower branch of (w), 
are well reproduced by the quasiparticle density of states. 
Above the Fermi level there is a peak in the quasiparti- 
cle density of states similar to that in the full spectrum 
but somewhat more pronounced. Above the Fermi level 
and below the upper peak there is a pseudo-gap region. In 
the free quasiparticle spectrum it is a definite gap. In the 
spectrum calculated from the direct NRG evaluation it ap- 
pears as a pseudo-gap, with rather small spectral weight 
just above the Fermi level. From the direct DMFT-NRG 
calculations, due to the broadening features introduced to 
obtain a continuous spectrum, it is not always possible 
to say definitively whether there is a true gap above the 
Fermi level or not. To resolve this question we can ap- 
peal to the renormalised perturbation theory to look at 
the corrections to the quasiparticle density of states aris- 
ing from the quasiparticle interactions. A calculation of 
the imaginary part of the renormalised self-energy Sai'-^) 
to order should be sufficient to settle this issue. The 
imaginary part of the second order diagram for the renor- 
malised self-energy in the limit T for > is given 
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Fig. 9. (Color online) The free local quasiparticle spectrum 
(dashed line) in comparison with DMFT-NRG spectrum for 
X = 0.9 and i7 = 3 for the spin-up electrons (upper panel) and 
spin-down electrons (lower panel). 



by 

CJ —UJ 

Imi:^2) (cj) ^irU'^Jdei Jde2 Po^<t(£i)/5o,-<t(w - £i + £2) 

(31) 





X P0-a{£2)d{t^ ~ £1 + £2), 



where po,o-(£) is the free quasiparticle density of states. 
The integration area is a triangle in the (ei, e2)-pla'iie as 
shown in figure [TOl 

To analyze the behavior oflmE^^ (w) in the regime |/io.T I < 
Lo < l/io.il we have to study where the integrand is non- 
zero taking into account that po^a-{£) = for l/iojl < 
£ < lAojil- The only non-zero contribution comes from 
the small shaded region in figure [101 which leads to the 
estimate, 

Im^(2)(^) ~ 7rU^poAO)po..-A-^)Po,-Mf^lr (32) 

When /io,T is small, which occurs when the lower edge of 
the gap in the quasiparticle density of states is very near 
the Fermi level, this contribution to the imaginary part 




Fig. 10. (Color online) Integration region in the (ei, £2)-plane 
for the imaginary part of the self-energy. The original triangle 
region (0, u, —u>) for integration in equation (|31|l is reduced in 
the gap region, I/^o.tI < ^ < lAoali to the small shaded region 
shown in the figure. 



of the renormalized self-energy will be finite but small. It 
decreases with u due the behavior of po,-it(— w). Based on 
this argument we conclude that there is a small, but finite 
imaginary part of the self-energy in the free quasiparticle 
gap 2Z\/2, when it lies above the Fermi level, giving rise 
to a finite spectral weight there. However, this spectral 
weight is very small close to the lower edge of the free 
quasiparticle density of states, when this edge lies only 
just above the Fermi level. 

In the lower panel of figure [S] we compare the Pi{uj) 
with Z|po.i(^)- We see that in this case also the quasi- 
particle density of states reproduces well the spectrum in 
the region of the Fermi level and the peak structure in 
the lower band, which is non-singular in this case. The 
position of the peak above the Fermi level is also well re- 
produced, but the peak in the free quasiparticle density 
of states is singular, whereas that in the DMFT-NRG re- 
sults is not. We would expect to lose this singularity in the 
free quasiparticle density of states once the quasiparticle 
scattering is taken into account and the renormalized self- 
energy is included. It is possible also that the peak above 
the Fermi level in the DMFT-NRG spectrum should be 
sharper, as there is some tendency for the broadening in- 
troduced in this approach to flatten peaked features in 
regions away from the Fermi level. The spectral weight 
in the pseudo-gap is even smaller than in the case for the 
spin-up electrons, particularly in the region of the gap that 
lies closest to the Fermi level. This is qualitatively in line 
with the conclusions based on the renormalized pertur- 
bation theory estimate of the effects of the quasiparticle 
scattering. 

We see very similar features in the spectra for the case 
U = 6 and also 10% doping shown in figure [TTl 
Here, the peaks near the Fermi level are a bit sharper. The 
observations made on the comparison of the quasiparticle 
and DMFT-NRG spectra apply equally well to this case. 
In addition to the low energy features charge peaks cor- 
responding to the Hubbard bands appear. The lower one 
can be identified in the full spectra, whereas the upper 
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Fig. 11. (Color online) The free quasiparticle spectrum 
(dashed line) in comparison with DMFT-NRG spectrum for 
X = 0.9 and U = 6 for the spin-up electrons (upper panel) and 
spin-down electrons (lower panel). 



Hubbard peak is not seen on the energy scale shown. The 
quasiparticle density of states does not contain informa- 
tion about these features at higher energy. 



4.2 fc-resolved Spectra 

We can learn more about the low energy single parti- 
cle excitations by looking at the spectral density of the 
Green's function G^^^{uj) in equation ([7]) for a given wave- 
vector k. With the self-energies Sai^) calculated within 
the DMFT-NRG approach all elements of this matrix can 
be evaluated. The local spectra and self-energies are spin- 
dependent in the doped broken symmetry state, however, 
the free quasiparticle bands [equation (fT9|) ] do not 

depend on the spin. Here, we focus on the diagonal part 
of Gfc ^ (w) corresponding to the A sublattice, 



Gfe,o-(w) 



(33) 



The weights of the quasiparticle excitations in this case 
depend on the spin corresponding to the sublattice prop- 



erties. We note that one can also analyze the quasiparticle 
bands differently, for instance, from the fc-resolved spectra 
and the diagonal form of Gj^ ai^)- The form of the quasi- 
particle bands remains unchanged then, but the weights 
differ and do not depend on the spin a in that case. 

We first of all look at the Fermi surface which is the 
locus of the fc-points at the Fermi level (w = 0) where 
the Green's function has poles. The conduction electron 
energy Sk-p at these point is given by 



(/iT-i:T(0))(A^i-^i(0)). 



(34) 



By Luttinger's theorem, the volume of the Fermi sur- 
face for the interacting system must equal that for the 
non-interacting system with the same density. As the self- 
energy depends only on oj, the two Fermi surfaces must 
also have the same shape, and therefore must be identical. 
The Fermi surface of the non-interacting system is given 
by Efcp = /io, where /ip is the chemical potential of the 
non-interacting system in the absence of any applied field 
for the given density. For this to be identical with that 
given in equation ()34p . 



(35) 



(/iT-i:T(0))(A^i-A(0))=A'^. 



We can check that this relation indeed holds from our 
results for Sai^^) and /ig-, independent of the value of U, 
or in the case of an applied staggered field, independent of 
the field value. This relation implies that the total number 
of electrons per site n can be calculated from an integral 
over the non-interacting density of states. 



po{uj)duj, 



(36) 



where in the hole doped case /io = — ^/a'tTh and fia- = 

To relate this result to the quasiparticle picture, we 
expand the self-energy in equation ()33|) to first order in 
u!, but retain the remainder term, 2J^{uj) as in equation 
P7p . The Green's function can be rewritten in the form. 



G'fc,<T(w) = J 



(37) 



where Co-(^) = + jj-o.a- — ^a{^)- We define a quasipar- 
ticle Green's function Gk,a('^) via ZcrGk,a{'^) ~ Gk.a-i^^)- 
The renormalized self-energy vanishes, Sai^^) — 0, for the 
free quasiparticle Green's function G^k ai'-^)' which can be 
separated into two independent branches of free quasipar- 
ticles, 

^(0),, ,^_ <(efc) , WLiek) 



- - El 



- - El 



(38) 



where ^ was defined in equation pop and the weights 
are given by 



Afl 



"3" 



(39) 
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This is similar in form to mean field theory, which would 
correspond to putting Zg- — I, and Ajj — Uirimi, where 
TTinif is the mean field sublattice magnetization. The spin 
dependent contribution in (|39p which arises from the sec- 
ond term is most marked in the region near the Fermi 
level. It should be noted that the quasiparticle excitations 
E'^ ^ and weights ^^-(fffc) here are defined by expanding 
the self-energy at w = 0. This is so that they correspond 
to the free quasiparticles in the renormalized perturbation 
theory which have an infinite lifetime. 

The spectral density p^,*'-' (to) for this free quasiparticle 
Green's function is a set of delta- functions, 

Pkii^) = - El^)+u^_{ek)6{cu - El_). (40) 

On the Fermi surface E^ _ = 0, which is consistent with 
the result for the Fermi surface given in equation (p4|) . 
Summing over k gives the local quasiparticle density of 
states in equation (^0)1 . We define the quasiparticle num- 
ber n as the integral of the sum of the spin up and spin 
down quasiparticle density of states up to the Fermi level, 



duj{LU + fl) (^/{iu + W^^A^ 



If we change the variable of integration to lo' , where 



(41) 



r,2 



^V^T^i V(^ + A)^ - 



the integration can be shown to be identical with that 
in equation ((36|) . using the fact that fiQ — —^/p-ip-i- We 
then have an alternative statement of Luttinger's theorem 
in the form n — n. This can also be found by summing 
both spin components in ()40|) . integrating over uj and then 
converting the fc-summation to an integral over the free 
electron density of states Pq{(-u). We can check in our nu- 
merical results that the relation in this form holds. The 
occupation number n can be calculated both from a direct 
evaluation of the number operator in the ground state, and 
also by integrating the sum of the spectral densities Pa{^) 
of the full local Green's function to the Fermi level. The 
value of n is similarly determined from the integral over 
the total quasiparticle density of states, Pai^)- AH three 
results were found to be in good agreement, to within one 
or two percent deviation at the most. 

Before discussing the fc-resolved spectra in detail we 
would like to ask what the spectral weight Wqp of a quasi- 
particle excitation at the Fermi level in the lower band is, 
such that the Green's function reads there 



(42) 



To calculate Wqp, we can not focus on the spin depen- 
dent local sublattice quantities, but have to sum over both 
sublattices or equivalently the two spin components. The 
reason for this is that the antiferromagnetically ordered 
state does not possess any net magnetization and has on 
average as many spin up polarized as spin down electrons. 



The division in the A and B sublattices is convenient for 
the DMFT calculations but somewhat artificial. In our 
case with hole doping the Fermi level lies within the lower 
band, which for the free quasiparticles is denoted by E'^ _. 
The corresponding weight on the Fermi surface defined by 
((34)) is then given by 



Wa 



where the average of the renormalized chemical potential 
p and the difference Ap were defined below equation (fT9|) . 
From the definition of Ap we can see that the second term 
in (j43p is spin rotation invariant. The spectral quasipar- 
ticle weight Wqp on the Fermi surface depends not only 
on the renormalization factors but also on the renor- 
malized chemical potentials po,a-- The same result for the 
weight (|43|) can be obtained from the diagonal form of 
Gfcg.(ti') and the spectral weight of the lower band. The 
weight Wqp corresponds to the spectral weight Z at the 
Fermi level as for example given in references 



The first term of the result for Wqp is like the arithmetic 
average of Za-. From figures [5] and [6] we can see that z^ > z^ 
and from figure [7] that /io,4 < /io.f < 0. Therefore the sec- 
ond term in (|43[) gives a positive contribution to the spec- 
tral weight. At the end of the section in figure flSl we show 
values of Wqp in comparison with the arithmetic average 
of z„. 

In order to understand better the properties of the 
quasiparticle bands, we now compare the quasiparticle 
spectrum with the fc-resolved spectral density pk^a-{uj) de- 
rived from the DMFT-NRG results. In figure [12] we make 
a comparison for the case of 12.5% doping with U = 3 
for the Green's function Gk.ai^) given in equation (I33p . 
Pk,cj{^) = —^'caGk.a{^^)l'^i where lo^ = oj + irj, with rj — > 
0, with that derived for the free quasiparticles, Zap'^\{uj) 
from equation The delta- functions of the free quasi- 
particle results are indicated by arrows with the height of 
the arrow indicating the value of the corresponding spec- 
tral weight. The plots as a function of oj are shown for a 
sequence values of Sk and, where the peaks in pk.a{^) 
get very narrow and high in the vicinity of the Fermi 
level, they have been truncated. It can be seen that the 
free quasiparticle results give a reasonable picture of the 
form of pfc^(j(w), particularly in the immediate region of 
the Fermi level. There is considerable variation along the 
curves in the way the overall spectral weight is distributed 
between the excitations below and above the pseudo-gap 
as a function of Sk- This is most marked in the region near 
the Fermi level for the spin-up electrons (upper panel) 
where most of the spectral weight is in the lower band and 
it is much reduced in the upper band, whereas the opposite 
is the case for the spin-down electrons. This is reflected 
in the expression of the quasiparticle weights u'^{ek) in 

equation ([5^ . For instance, uL(efc) corresponding to the 
lower band E'^ _ becomes maximal near the Fermi en- 
ergy, whereas u\{ek) goes to zero there. The finite width 
of the quasiparticle peaks in pk^ai^) can be described by 
a RPT, when we take into account the renormalized self- 
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Fig. 12. (Color online) The spectral density ph,a{i^) for the 
spin-up electrons (upper panel) and spin-down (lower panel) 
plotted as a function of uj and a sequence of values of Sh for 
[7 = 3 and 12.5% doping. Also shown with arrows are the 
positions of the free quasiparticle excitations, with the height 
of the arrow indicating the corresponding weight. 



energy S„{lo) in equation (|37p . If we, for instance, use 
the the second order approximation in C/, which was illus- 
trated in the last section ((3T|) . we get a similar behavior 
for small to as seen for pk^^iuj) in figure fT2l 

From the positions of the peaks in the Pk.<j{<^) spectra 
we can deduce two branches of an effective dispersion Ek,± 
for single particle excitations and compare it with the ones 
for the free quasiparticles ^ . We give the results for 
J7 = 3 in figure [131 

It can be seen that i?^ _ tracks the peak in the lower 
band closely over a wide range of Sk, —1.5 < Sk < 1-5 
(note the bandwidth W = A). This is not the case in the 
upper band, where _,_ tracks the peak closely only in 
the lowest section that lies closest to the Fermi level. As 
one can see from the dotted line the Fermi level lies in 
the lower band and intersects the lower band twice. This 
corresponds to the two values with opposite sign as 
can be see from equation 

The corresponding results for the fc-rcsolved spectra 
for J7 = 6 and also 12.5% doping are shown in figure [TH 
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Fig. 13. (Color online) A plot of the peak positions Ek,± in 
the spectral density pk^ai^) (full line) as a function of Sk for 
U = 3 and 12.5% doping compared with the free quasiparticle 
dispersion (dashed line). 




Fig. 14. (Color online) The spectral density pk^„{uj) for the 
spin-up electrons (upper panel) and spin-down (lower panel) 
plotted as a function of u} and a sequence of values of Sk for 
U — 6 and 12.5% doping. Also shown with arrows are the 
positions of the free quasiparticle excitations, with the height 
of the arrow indicating the corresponding weight. 
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In order to compare well with the case C/ = 3 we have 
chosen an identical range for oj and £k , although the large 
spectral peaks near the energy are very close together in 
this presentation. It can be seen that the overall features 
are very similar to those seen for U — i. For the spin 
up spectrum (upper panel) the peaks for the lower band 
have most of the weight near the Fermi energy, whereas 
the upper band is suppressed there, and vice versa for 
the opposite spin direction. The lower bands are tracked 
well by the free quasiparticles, and we can see that the 
bands for the larger value of U are significantly flatter. 
This is also clearly visible in the following figure [151 where 
we again compare the quasiparticle band with the peak 
position of the full spectra. On the range shown the lower 
band Ef. _ completely coincides with the free quasiparticle 
band Ei_. 
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Fig. 15. (Color online) A plot of the peak positions Eh.± in the 
spectral density pk.aiuj) (full line) as a function of Sh for U = Q 
and 12.5% doping compared with the free quasiparticle disper- 
sion (dashed line). On the range shown the lower band Ek,- 
completely coincides with the free quasiparticle band E^ _ . 



From the fc-resolved spectra in figures [T^] and [H] we can 
also extract the width of the quasiparticle peak Z\qp in 
the spectral density Pk.a-i^) (majority spin a =t). Its in- 
verse 1/^qp gives a measure of the quasiparticle lifetime. 
The results for Z\qp for the lower band £'fc._ for the two 
cases [/ = 3, 6 and 12.5% doping are shown in figure [TBI as 
function of Cfc. 

This plot brings out more clearly the feature that can 
be seen already in figures [T2] and [14] (upper panel) that 
the width increases sharply when we move away from the 
Fermi level and the values for the width Z\qp for U = 6 
are significantly smaller than those for U — S. This is in 
line with the fact that the local quasiparticle interaction 
U is smaller for the larger value of the bare interaction U 
as commented on earlier. The free quasiparticle picture is 
therefore even more appropriate in the case with stronger 
interaction. To numerical accuracy the width vanishes at 



and is finite for the interval Si 



< Ek < £k which lies 
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Fig. 16. (Color online) A plot of the width of the peaks Z\qp in 
the spectral density of the majority spin pfc.|(a;) as a function 
of Eh for ?7 = 3 (dashed line) and U = 6 (full line) and 12.5% 
doping. 



Another quasiparticle property, the effective mass en- 
hancement m* /m, can be extracted by calculating the 
derivative of _ in ([19]) with respect to Efc, which yields 
when evaluated at the Fermi energy ([M]) . 



1 



(44) 



within the lower band but above the Fermi level. 



The effective mass enhancement therefore does not only 
depend on z^, but also on the renormalized chemical po- 
tentials /io,(T. The general trend for m*/m as function of 
U can be seen in figure [T71 for the case of 7.5% doping. 
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Fig. 17. The ratio m*/m according to (|44p plotted over a 
range of /U for 7.5% doping. 



The effective mass increases sharply for large U as the 
hole motion is energetically more costly in the ordered 
background. The fact that the lower band for [/ = 6 seen 
in figure [15] is flatter than in the case U = 3 m figure [13] 
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can be attributed to the larger effective mass. We find a 
similar behavior of m* /m as function of U for different 
filling factors from the ones shown in figure 1161 The trend 
is that the effective mass enhancement is less pronounced 
for larger doping, which is intuitively understandable by 
the quasiparticle motion in an ordered background. 

In the DMFT framework for the paramagnetic state 
as well as the case with homogeneous magnetic field, the 
quasiparticle spectral weight Wqp and the inverse of the ef- 
fective mass enhancement m/m* can be described simply 
by the renormalization factor z„. In figure [18] we show a 
comparison of the spectral quasiparticle weight Wqp (|43p . 
the arithmetic, {z-\ + zx)/2, and geometric, y^zfzj", aver- 
age of the renormalization factors, and the inverse of the 
effective mass, m/m*, (|44)) for [/ = 3 for various dopings. 



0.7 




0.05 0.1 0.15 

5 

Fig. 18. (Color online) Comparison of the spectral quasiparti- 
cle weight lUqp from equation (|43|l . the arithmetic, (z| -(-2;|)/2, 
and geometric, ^z-^zi, average of the renormalization factors, 
and the inverse of the effective mass, m/rn' from equation (|44p . 
for (7 = 3 and a range of dopings. 



As seen in this case with antiferromagnetic symmetry break- 
ing these quantities take a different form (|43)l and (|44p 
and have distinct values. For different values of U the 
behaviour is qualitatively similar. As a first approxima- 
tion the quasiparticle spectral weight Wqp corresponds to 
the arithmetic average of the renormalization factors Za, 
whilst m/m* relates to the geometric average. In general, 
one can, however, not omit the dependence on the renor- 
malized chemical potential as it gives a significant contri- 
bution as can be seen in figure [TSl This can be understood 
for example for the limit of zero doping. The system then 
becomes an antiferromagnetically ordered insulator with 
spectral gap. The weights z^, tend to finite values, but the 
effective mass must diverge. This is found in equation (j44p 
since /io.T ~* for (5^0, and the trend can be seen in 
figure [m 



5 Conclusions 



We have studied the field induced and spontaneous anti- 
ferromagnetic ordering in the hole doped Hubbard model 
with DMFT-NRG calculations at T = 0. A phase diagram 
separating antiferromagnetic and paramagnetic solutions 
for different values of doping and interactions U ranging 
from zero to about 1.5 times the bandwidth W has been 
established and is in agreement with earlier results by 
Zitzler et al. [lOj . Our main objective has been to ana- 
lyze the properties of the quasiparticle excitations in the 
metallic antiferromagnetic state. We presented two differ- 
ent ways of calculating the parameters z„ and /io,(j , which 
define the renormalized quasiparticles, and the two sets 
of results have been shown to be in agreement. We have 
also been able to deduce the effective on-site quasiparticle 
interaction U from the NRG low lying excitations. The 
low energy properties of the local spectral function can 
be understood in terms of the free quasiparticle picture. 
We have used the second order perturbation expansion in 
powers of U to estimate the spectral weight in the pseudo- 
gap region above the Fermi level. 

We have been able to compare the position of the 
peaks found in the fc-dependent spectral functions with 
the dispersion relation for the free quasiparticles. The free 
quasiparticle dispersion gives a very good fit to the posi- 
tion of these peaks in the lower band which intersects the 
Fermi level. The quasiparticle lifetime, as deduced from 
the widths of the peaks in the spectrum, increases for 
stronger interactions. This is consistent with the fact that 
the on-site quasiparticle interaction U , which gives the 
quasiparticles a finite lifetime, decreases with increase of 
U in the same range. We have also shown how the spec- 
tral quasiparticle weight at the Fermi level Wqp and the 
effective mass can be deduced from the parameters z^ and 
/io,(T- The effective mass is found to increase with the in- 
teraction, and it diverges in the limit of zero doping whilst 
Wqp remains finite. 

We have found that Luttinger's theorem for the total 
electron density in the antiferromagnetically ordered state 
holds within the numerical accuracy for the range of dop- 
ings and interactions studied. This is a further indication 
that many aspects of Fermi liquid description may hold in 
situations with symmetry breaking. 

It is not easy to make a direct comparison of our re- 
sults with earlier work [3] analyzing the quasiparticle exci- 
tations in an metallic antiferromagnet as these have been 
mainly based on the t — J-model for one or two holes 
in a finite cluster. However, at a semiquantitative level, 
the overall trend in our results seems to be similar to the 
results surveyed by Dagotto, where the effective quasipar- 
ticle bandwidth Woff is found to decrease with decreasing 
J. This is line with our results if we identify W^s ~ m/m* 
and J ^ /U (see figure [17)) . Our values for the spectral 
quasiparticle weight Wqp are qualitatively similar to those 
presented as the wavefunction renormalization Z in the 
review article by Dagotto (see fig 27 [5]), and also the 
ones reported more recently [32] • 
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